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Abstract 

A package for pattern recognition and PID by COMPASS RICH-1 has been developed and used for the analysis of 
COMPASS data collected in the years 2002 to 2004, and 2006-2007 with the upgraded RICH-1 photon detectors. It 
has allowed the full characterization of the detector in the starting version and in the upgraded one, as well as the 
PID for physics results. We report about the package structure and algorithms, and the detector characterization 
and PID results. 
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1. Introduction 

COMPASS PQ is a fixed target experiment at 
CERN SPS. The COMPASS physics program is 
focused mainly on the study of the nucleon spin 
structure and on hadron spectroscopy. COMPASS 
has collected data in 2002-2004 and 2006-2007 
using a muon beam of 160 GeV/c. The appara- 
tus includes a polarized target for spin structure 
studies, a large acceptance double spectrometer, 
electromagnetic and hadron calorimeters and muon 
filters [2]. Charged hadron identification (PID) is 
performed by means of a focusing Ring Imaging 
CHerenkov detector of large dimensions, the COM- 
PASS RICH-1 0. RICH-1 is a Cherenkov imaging 
detector with large angular acceptance (±180 mrad 
vertical, ±250 mrad horizontal). The radiator gas 
is C4F10 at atmospheric pressure and the typical 
particle path length in the radiator is 3 m. Im- 
age focusing is obtained by 2 spherical surfaces, 
formed by hexagonal and pentagonal mirror ele- 
ments, resulting in a reflecting surface with a total 
area larger than 20 m 2 . Images are collected on 
two sets of photon detectors, placed above and be- 
low the detector acceptance region. Until 2004, the 
photon detectors used were eight MWPCs with seg- 
mented Csl photo-cathodes, covering a surface of 
5.2 m 2 . In 2006, an important upgrade of the pho- 
ton detectors system was completed. The upgrade 
is based on two complementary techniques for the 
central and peripheral regions of the photon de- 
tectors. The central part, corresponding to four of 
the 16 MWPC photo-cathodes, has been equipped 
with multi-anode photomultiplicr (MAPMT) with 
UV extended glass window (Hamamatsu R7600- 
03-M16) coupled to individual telescopes of fused 
silica lenses and a fast digital read-out system [4]; 
in the peripheral part, a new readout system [5], 
based on the APV chip [6], almost dead-time free 
and characterized by improved time resolution, has 
been installed on the already existing Csl MWPCs. 

From the point of view of the software data han- 
dling, the upgrade results in two main consequences: 
- The geometry and the spatial resolution of the 
central and peripheral detectors are different: the 
size of the pads in the peripheral part is 8 x 8 mm 2 , 
while in the central part the pseudopads (the 
MAPMT pixels projected through the optical 
telescope onto the plane previously housing the 
Csl photo-cathodes) have a dimension of about 
12 x 12 mm 2 ; the different spatial resolution is one 



of the elements that requires different detector 
characterization in the two parts. 
- The spectrum of the detected Cherenkov photons 
is different in the central part (200-750 nm) re- 
spect to the spectrum of the peripheral part (165- 
200 nm); this implies that the average effective 
value of the radiator refractive index is different 
according to the type of the photon detector in- 
volved; the mean values of the effective refractive 
index, as evaluated offline from a data sample, are 
respectively: n - 1 ~ 0.001345 for the MAPMT 
detectors andn-1 - 0.001528 for the Csl MWPC 
detectors. The difference in the effective refractive 
indexes implies also that the rings detected in the 
MAPMT and the Csl region, have different radii: 
the radius of the rings detected in the peripheral 
region are about 30% larger than those detected 
in the MAPMT region. 

2. RICH-1 software 

The analysis of the RICH data can be divided 
in two main steps: the ring reconstruction and the 
particle identification. The aim of the RICH in the 
experiment is to give an answer to the latter point, 
nevertheless a complete understanding of the RICH 
performances is needed before any attempt to iden- 
tify particles, and for this purpose the ring recon- 
struction is largely used. 

The reconstruction of the RICH-1 data is per- 
formed with RICHONE, a dedicated class of the 
COMPASS reconstruction programme, CORAL. As 
starting point of the reconstruction, the RICH dig- 
its in both the MAPMT and the Csl MWPC part 
are used; the hits are rearranged in clusters in the 
Csl MWPC part only, since the MAPMT crosstalk 
is completely negligible. Using the cluster position 
on the detector surface, and assuming as emission 
point the middle point of the particle trajectory in- 
side the radiator, the photon emission angles respect 
to the particle trajectory (Ochtfich) are computed 
using the Hough transformation, following a recipe 
well known in the literature [7J. The trajectory is 
provided by the tracking package in the CORAL re- 
construction programme. 

The ring reconstruction is mainly used for the 
analysis of the RICH performances. For all the clus- 
ters contained in a region of 70 mrad around the 
particle trajectory, the emission angles (Ochj^Ch) 
are computed. Since a ring is characterized by a 
fixed value of 9c hi a peak in the 9cn distribution is 



2 



searched for through a scan with a fixed size window 
of ±3(T p h, where a p h is the single photon resolution. 
For the ring reconstruction, the Cherenkov angles 
of the photons detected in the MAPMT region are 
normalized to the angles in the Csl MWPC region, 
so as to combine, for rings not completely included 
in a single detector region, the information coming 
from the two different regions. One of the main ad- 
vantages of the simple recipe described above is the 
automatic handling of events of different complex- 
ity, as is the case of the split rings. In fact, due to 
the RICH-1 detector architecture, photons emitted 
by particles scattered at small vertical angles, illu- 
minate both mirror surfaces, and are then reflected 
partially onto the upper and partially onto the lower 
photon detector set, forming a full ring image in a 
detector set and a partial ring in the other one. 

The particle identification algorithm relies on a 
likelihood function built from all the photons associ- 
ated to the particle in the fiducial region of 70 mrad. 
A Poissonian probability is used to describe the ex- 
pected number of detected photons (both signal and 
background contribution) as a function of the pa- 
rameters of the particle trajectory. The likelihood is 
computed for five mass hypothesis (e, v, tt, K,p) and 
for the background, corresponding to the hypothe- 
sis of absence of signal. For a given mass M, and for 
N photons, the likelihood expression reads: 
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Fig. 1. Number of photons per ring as a function of 
the Cherenkov angle, for rings detected in the MAPMT 
part only; the curve is a fit with a function of the type 
No sin 2 (Ac/,). The number of photons at saturation after 
background subtraction is around 56. 
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Fig. 2. Standard deviation of the 8 r i ng — 8-n distribution, 
for particles identified as pions, as a function of the particle 
momentum; ring detected in the MAPMT detector only. 
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The background term, b(9j, <pj), is evaluated pad 
by pad from a map of the integrated cluster distri- 
bution, taken from the data themselves. Sm and B 
are the signal and background photons integrated 
on the particle fiducial region, defined above. The 
signal term sm(9j-, <pj) is taken as a Gaussian, with 
mean value the Cherenkov angle evaluated from the 
kinematics corresponding to the mass M and with 
the width of the single photon resolution: 
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Sq is the expected number of photons evaluated us- 
ing the Frank and Tamm equation; the e term is the 
probability of detecting a photon, taking into ac- 
count the dead zones in the detector. 



3. Detector characterization 

RICH-1 has been fully characterized in both the 
old and the upgraded version, monitoring in opera- 
tive conditions all the relevant observables: the num- 
ber of detected photons, the angular resolution and 
the PID efficiency. 

The number of signal photons emitted at satura- 
tion has been evaluated through a fit of the number 
of photons per ring as a function of the Cherenkov 
angle, using a function of the type N sin 2 (9ch)- 
The number of detected photons at saturation is 
around 14 before the RICH-1 upgrade and, after the 
upgrade, in the peripheral regions, while it has in- 
creased up to 56 (fig. [J) in the upgraded RICH, cen- 
tral regions. Moreover, the ratio of the signal to the 
background in each ring is increased because the 
number of background clusters per ring has been re- 
duced by a factor ^30% in the central part. 

The single photon resolution is evaluated from the 
width of the distribution 6 p h — where 9 p h is the 
Cherenkov angle of a photon belonging to the ring 
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Fig. 3. PID efficiency for K — , evaluated on a sample of 
particles from <j> decay, as a function of the particle polar 
angle. The two sets of data points correspond to the old and 
to the upgraded versions of RICH-1. 
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Fig. 4. Same as fig. [3] for K+ 

and 9 n is the angle for the pion mass hypothesis. In 
RICH-1 before the upgrade, the mean value of the 
resolution is around 1.2 mrad, while the resolution 
on the ring angle is around 0.6 mrad. The two num- 
bers do not scale with the square root of the number 
of detected photons per ring due to the large back- 
ground contribution in each reconstructed ring, di- 
luting the signal. After the upgrade, the single pho- 
ton resolution in the central part is around 2 mrad, 
while the ring angular resolution less than 0.3 mrad 
(fig. [2]); the better scaling of the ring resolution with 
the number of photons is due to the fact that the 
background contribution per ring is small. 

The PID efficiency has been evaluated selecting 
with kinematics criteria a sample of exclusively pro- 
duced <p mesons; the sample has a purity of about 
90%. The efficiency has been evaluated separately 
for K + and K~ , using the positive or negative track 



coming from the <fi decay. In fig. [3] and [H the K ef- 
ficiency is shown as a function of the particle po- 
lar angle, both for the old and for the upgraded 
RICH. The curve corresponding to the old RICH 
shows clearly that the efficiency increases at large 
polar angle, since at small angles it is limited by the 
presence of an important background coming from 
the muon beam halo. The impact of the upgrade on 
the RICH efficiency is clearly visible from the corre- 
sponding efficiency curve: also at small polar angle 
the efficiency is above 90%. 

4. Conclusions 

The RICH-1 reconstruction software is based on 
a simple recipe, allowing an easy handling of events 
of different complexity. The ring pattern recogni- 
tion and the PID are fully independent and are used 
for different purposes: detector characterization and 
data analysis, respectively. The complete RICH-1 
characterization, both in the old and in the upgraded 
version, has been presented and confirms a large 
improvement of the RICH-1 performance in its up- 
graded version. 
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